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The idea of partons and the utihty of using hght-cone gauge in QCD are 
introduced. Saturation of quark and gluon distributions are discussed using 
simple models and in a more general context. The Golec-Biernat Wiisthoff 
^ , model and some simple phenomenology are described. A simple, but real- 

Tj- I istic, equation for unitarity, the Kovchegov equation, is discussed, and an 

^ I elementary derivation of the JIMWLK equation is given. 

1 Introduction 

', These lectures are meant to be an introduction, and an overview, of par- 

I I ton saturation in QCD. Parton saturation is the idea that the occupation 

numbers of small-x quarks and gluons cannot become arbitrarily large in the 
^ ■ light-cone wavefunction of a hadron or nucleus. Parton saturation is an idea 

which is becoming well established theoretically and has important applica- 
^ ! tions in small-x physics in high-energy lepton-hadron collisions and in the 

5j I early stages of high-energy heavy ion collisions. The current experimental 

situation is unclear. Although saturation based models have had considerable 
phenomenological success in explaining data at HERA and at RHIC more 
complete and decisive tests are necessary before it can be concluded that par- 
ton saturation has been seen. These lectures begin very simply by describing, 
through an example, some important features of light-cone perturbation the- 
ory in QCD, and they end by describing some rather sophisticated equations 
which govern light-cone wavefunctions when parton densities are very large. 



^Lectures given at the Cargese Summer School, August 6-18, 2001 
^This work is supported in part by the Department of Energy 
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FIG.l 

2 States in QCD Perturbation Theory 

In light-cone perturbation theory states of QCD are described in terms of 
the numbers and distributions in momentum of quarks and gluons. The 

essential features of light-cone perturbation theory that will be needed in 
these lectures can be illustrated by considering the wavefunction of a quark 
through lowest order in the QCD coupling g. One can write 

\i;,>=N\p> + Y, E d'krx{k)\p-k;k{X,c)> (1) 

A=± c=l 

where |p > is a free quark state with momentum p, I'^p > is a dressed quark 
state, and \p — k; k{X, c) > is a state of a quark, of momentum p — k, and 
a gluon of momentum k, helicity A and color c. We suppress quark color 
indices which will appear in matrix form in what follows. Wc label states 
by momenta p+ = ^ {po + P3),Pi,P2, and d^k = dk^d^k = dk^dkidk2. 
Recall that in light-cone quantization momenta P+ and P — (Pi,P2) are 
kinematic while P_ plays the role of a Hamiltonian and generates evolution 
in the "time" variable x+ — ^ (xq + Xs). For an on-shell zero mass particle 
P- — p^/2p+. in (1) is a normahzation factor. Eq.(l) is illustrated in Fig.l. 



■0^ is determined from light-cone perturbation theory to be 
Sip-p)Uk)^ + 

with 

Hi^gJ d'xq{xh,{^)q{x)A;{x) (3) 
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where the gluon field is 



Alix) f [6;(A;)a^(fe)e^^-^-^^+— + h.c] (4) 

with 

[ai{k'),a^^{k)]^5^.x5e,5\k' -k). (5) 

In (2) and (5) 5^{p) = S{p+)S'^{p) while d^x = dx^cPx in (3). It is useful to 
imagine the calculation being done in a frame where p+ is large and p = in 
which case 

In the soft gluon approximation k^/p^ << 1 and thus k_ >> {p—k)_ so that 
only /c_ need be kept in the denominator in (2). In addition, in light-cone 
gauge, = 0, the polarization vectors can be written as 

6i(fc) = (6^,6^6^) = (0,^,6^) (7) 

and, because of the 1/k^ term, only need be kept in (2) in the soft gluon 
approximation. Using (3)- (7) in (2) one finds 

- i^)29%^ I, \ , (8) 
2 k y'(27r)32A;+ 

Problem 1(E): Using the formula U{p — k)^fJJ{p) ~ 2p^gn_ for high mo- 
mentum Dirac spinors derive (8). 

3 Partons 

Define the gluon distribution of a state \S{p) > by 
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xGs{x,Q^) = E/ d^kx5{x-k+/p+)e{Q^-k^) < S{p)\a'^{k)al{k)\S{p) > . 

A,c 

(9) 

The meaning of xGs is clear. xGs{x, Q'^)dx is the number of gluons, having 
longitudinal momentum between xp+ and + localized in transverse 

coordinate space to a region Ax±_ ~ l/Q, in the state |<S'(p) > . For a quark, 
at order g"^, one finds from (1) 

xG,{x, Q^)^T.j d^'kxSix - k+/p+)e{Q^ - k^)r^{k)rx{k). (10) 
Using (8) in (10) one finds 

xG,{x, Q') = Eyy/(^^^ ^^(^ - k^/p^MQ' - (11) 

Using Tic ^ ^ — Cp — ^2N^ introducing an infrared cutoff, for the 
transverse momentum integral in (11) gives 

xG^{x, g2) = ^£n(g2/^2). (12) 

TT 

We note that if xG{x, Q^) = 3xGq{x, Q^) is taken one obtains a result for the 
proton which is not unreasonable phenomenologically for x ~ 10~^ — lO"-*^ 
and moderate if n is taken to be 100 MeV. 



4 Classical Fields 

One can associate a classical field with gluons in the quark. 

Af\x) = / < ^I^MtixMp > . (13) 
Using (1) and (8) in (13) one finds 



^rv) - 1 (14) 



or 



with 



A^rKx) = / ^,e-'<''4"\k) (15) 
^"(1=) = y f . (16) 



In (14) the A;+ integration goes from — oo to +oo. The region > comes 
from \%l)p > consisting of a bare quark and a gluon of momentum k while the 
region A;+ < comes from \ipp> > consisting of a bare quark and a gluon. 

Problem 2(E): Take ^ = in (16) and show that 

Problem 3(E): Show that Gq, as given in (10), can also be written as 

xGq{x,Q')^ I -^^5{x_k^/p^)Q{Q' -k')Y.[Ar\k)f2K- 
i^^r i.c 



5 Why Light-Cone Gauge is Special 

In order to understand why light-cone gauge plays a special role in describing 
high-energy hadronic states a simple calculation of the near forward high- 
energy elastic amplitude for electron-electron scattering in QED is useful. 
The graph to be calculated is shown in Fig. 2 and we imagine the calculation 
being done in the center of mass frame with p being a right-mover (p+largc) 
and pi being a left-mover {pi_ large). First we shall do the calculation in 
covariant gauge and afterwards in light- cone gauge. 



In covariant gauge the photon propagator, the k-line in Fig. 2, is D^,^ 
^g^u and the dominant term comes from taking g^i, —>■ = 1 giving 



^ ^ . U{p-kh+U{p) . U{pi + kh^U{pi) -4 

2p+ 2pi_ ^ ^ 
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FIG.2 

Using U{p — A;)7_|_t/(p) pa 2p+ and U{pi + k)'^-U{pi) ^ 2pi_ 
one finds 

T = ^ (18) 

where we have used k"^ = 2A;+A;_ — ^ —h^ since both and A;_ are 
required to be small from the mass shell conditions {p — k)"^ = (pi + A;)^ = 
for the zero mass electrons. 

Now suppose we do the calculation in hght-cone gauge — where the 
propagator is 

j.-w-gfo^- ""^;"';"*'' ] (19) 

where r) • V — V+ ior any vector V^. Now the dominant term comes from 
taking D_i giving 

_ . U{p - k)^+U{p) . ^(pi + fc)7-fc^(pi) -i 
^"'^ 2p+ 2pi_ ^^^^ 

Using f7 (pi + A;)7 ■ ^ U{pi) = k^ one finds 

T=^^. (21) 

Now {pi + ky = gives 2pi_A:+ ^ so that (18) and (21) agree as expected. 

The result (18) comes about in a natural way ^'^p^'^J+^'^p^ and ^(Pi+^^^h-'^fa) 
are the classical currents, equal to 1, of particles moving along the light- 
cone while the 1/k^ factor is just the (instantaneous) potential between the 
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charges. However, when one uses = hght-cone gauge the dominant 
part of the current for left moving particles is forbidden and one must keep 
the small transverse current fc/2pi_. However, the smallness of the current is 
compensated by the factor 1/ A;+ in the light-cone gauge propagator. In coor- 
dinate space the l/k+ comes about from the potential acting over distances 
x„ ~ l/fc+ = so that the potential is very non-local and non-causal. 
By choice of the ie prescript ion 0] for l/fc+ one can put these non-causal 
interactions completely before the scattering (initial state) or completely af- 
ter the scattering (final state). In Problem 2 the potential Ai exhibits this 
non-causal behavior with our choice of ie placing the long time behavior in 
the initial state, the B(— x_) term. 

6 High Momentum Particles and 
Wilson Lines^"^ 

Suppose a quark of momentum p is a high-energy right mover, that is p+ >> 
P-,p. Then so long as one does not choose to work in A_ = light-cone 
gauge the dominant coupling of gauge fields to p are classical (eikonal) when 
p passes some QCD hadron or source. To be specific suppose p scatters on 
a hadron elastically and with a small momentum transfer. We may view the 
interactions as shown in Fig. 3 for a three-gluon exchange term. Of course 
to get the complete scattering one has to sum over all numbers of gluon ex- 
changes. Call S{p^,b) the S-matrix for scattering of the right moving quark 
on the target. (We imagine that the target hadron has large gluon fields 
making it necessary to find a formula which includes all gluon interactions.) 
Although the right moving quark has sufficient momentum so that it moves 
close to the light-cone we do assume that the momentum is not so large that 
higher gluonic components of the quark wavefunction need be considered. 
Since the probability of extra transverse gluons being present in the wave- 
function is in general proportional to ay, with y being the quark rapidity, we 
suppose ay « 1. 

The graph shown in Fig. 3 can be written as 
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where 



Mabc{ki,k2) = J d^xid^X2d^xe'(*^-*^i-*^^>^+^'=2-^2+*-^i 

• < f\TAl{x)A'4x2)Al{x,)\i > . (23) 
It is straightforward to evaluate the { }-term in (22) and one gets 

{ } = igT'^^gT''^gT^-^ \ (24) 

ki- + ie {ki + k2)- + t€ 

Now do the d'^kidki^ and d'^k2dk2+ integrals followed by d^Xidxi^ and d^X2dx2- 
This sets x_i — X2 — x and xi- — X2- — X-. The only non-zero term in the 
time-ordered product is proportional to ©(x_ — X2-)Q{x2- — xi^) and this 
factor along with the exponentials in Xi_ , X2- and a;_ allow the dki- and 
dk2- integrals to be done over the poles in (24). We get finally 

/oo rx-\. 
dx+A^{h,x+) / dxi+A_{h,X2+) 
-oo J —oo 

dxi+A^{h,xx+)\i> (25) 

-oo 
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where = T°'A'^^ and we have suppressed the variable x_ = in the A' s in 
(25). 

The general term is now apparent, and one has in the general case 

Sfi{p+A) =< /|Pe^^-/'-°°-''"+^-(-'"+)|z > (26) 

where P denote ordering of the matrices A where A's having larger 

values of x+ come to the left of those having smaller values. 

Problem 4(M): Show that the time orderings different from Q{x^—X2+)Q{x2+ — 
X14.) give no contribution to (22). 

Finally, a word of caution in using (26). In general there are singularities 
present when two values of x+, in adjoining A's, become equal, although 
in many simple models the x+— integrations are regular. If singularities 
in the x+ integrations arise it is generally possible to extract the leading 
logarithmic contributions to the scattering amplitude by carefully examining 
the singularities 0. A detailed discussion of this is, however, far beyond the 
scope of these lectures. 

7 Dual Descriptions of Deep Inelastic Scat- 
tering; Bjorken and Dipole Frames 

Particular insight into the dynamics of a process often occurs by choosing a 
particular frame and an appropriate gauge. Indeed the physical picture of a 
process may change dramatically in different frames and in different gauges. 
In a frame where the parton picture of a hadron is manifest saturation shows 
up as a limit on the occupation number of quarks and gluons, however, in a 
different (dual) frame saturation appears as the unitarity limit for scattering 
of a quark or of a gluon dipole at high-energy 

To see all this a bit more clearly consider inelastic lepton-nucleon scatter- 
ing as illustrated in Fig.4 where a lepton emits a virtual photon which then 
scatters on a nucleon. We suppose = —q^ is large. The two structure 
functions, Fi and F2, which describe the cross section can depend on the 
invariants and x = ^r^. To make the parton picture manifest we choose 
A+ = light-cone gauge along with the frame 
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and 

<1 = (qo,g,(lz = 0), 

and where P — > cxd. This is the Bjorken frame. We note that qo = ^ goes to 
zero as P — > oo so that the virtual photon momentum is mainly transverse 
to the nucleon direction. This last fact means that the virtual photon is a 

good analyzer of transverse structure since it is absorbed over a transverse 
distance Ax± ~ l/Q. Since Aa;^ is very small at large Q the virtual photon 
is absorbed by, and measures, individual quarks. 

Problem 5(H): Use the uncertainty principle to show that the time, At, 
over which j*{q) is absorbed by a quark is Ar ~ 2xP/Q'^. You may assume 
that ^ IQ^ << 1 where the process of absorption of the photon is illustrated 
in Fig. 5. 



The result of problem 5 and of the result that Axj. ~ l/Q motivates the 
formula 
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F2{x, Q') = Yl Q') + ^^/(^' Q')] (27) 

/ 

which says that the structure function F2 is proportional to the charge 
squared of the quark, having flavor /, absorbing the photon and propor- 
tional to the number of quarks having longitudinal momentum fraction x 
and localized in transverse coordinate space to a size 1/Q. F2 is given in 
terms of the longitudinal and transverse virtual photon cross sections on the 
proton as 



F2 = —. Q'[(Jt + (Jl]. (28) 

Eq.(27) is the QCD improved parton model. In more technical terms is a 
renormalization point which in the parton picture is a cutoff of the type given 
in (9) for gluon distributions and here occurring for quark distributions. 
Now consider the same process in the dipole frame pictured in Fig. 6 where 



q=i\/q'-Q',0,0,-q) 
11 




FIG.6 



and where q/Q » 1 but where q is fixed as x becomes small so that most 
of the energy in a small-x process, and it is only for small-x scattering that 
the dipole frame is useful, still is carried by the proton. Now we suppose 
that a gauge different from = is being used, for example a covariant 
gauge or the gauge = 0. In the dipole frame the process looks like 7* — *• 
quark — quark followed by the scattering of the quark- ant iquark dipole on 
the nucleon|0, |, ||. The splitting of 7* into the quark- ant iquark pair is 
given by lowest order perturbation theory while all the dynamics is in the 
dipole-nucleon scattering. In this frame the partonic structure of the nucleon 
is no longer manifest, and the virtual photon no longer acts as a probe of the 
nucleon. 

Equating these two pictures and fixing the transverse momentum of the 
leading quark or ant iquark one has0 

^d{xqf + xqf) r (Pxi(Px2 , ^ , f* i n\ 1 f i n\ 

■ e-'^<^^-^^'>[S\x,)S{x,) - S\x,) - Six,) + 1], (29) 
where the 7* wavefunction is 

(30) 
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We have written (29) for a fixed impact parameter, b, which is a (suppressed) 
variable in S. To make the identification exhibited in (29) requires that the 
struck quark shown in Fig. 5 not have final state interactions, and this requires 
the special choice of ie's in the light-cone gauge, = 0, used to identify the 
quark distributions. With this choice of ie's the lefthand side of (29) refers 
to the density of quarks in the nucleon wavefunction while the righthand side 
of (29) refers to the transverse momentum spectrum of quark jets produced. 
This identification is only possible when final state interactions are absent. 
Finally, it is not hard to see the origin of the final factor on the righthand side 
of (29). This is just [S''{x2) — l][S'(a;i) — 1], the product of the T— matrices 
for the dipole Xi in the amplitude and the dipole X2 in the complex conjugate 
amplitude. 

8 Quark Distributions in a large Nucleus; Quark 
Saturation at One-Loop 

In general it is very difficult to evaluate the S— matrices appearing in (29). 
However, if the target is a large nucleus one can define an interesting model, 
if not a realistic calculation for real nuclei, by limiting the dipole nucleon 
interaction to one and two gluon exchanges. That is, we suppose the dipole- 
nucleon interaction is weak. Despite this assumed weakness of interaction the 
dipole nucleus scattering can be very strong if the nucleus is large enough. 
We begin with the term S{xi) in (29), an elastic scattering of the dipole 
on the nucleus in the amplitude with no interaction at all in the complex 
conjugate amplitude. This term is illustrated in Fig. 7. 



Now |S'(xi)p is the probability that the dipole does not have an inelastic 
interaction as it passes through the nucleus. We can write 

|5(xOr = e-^A (31) 

where L = 2\J E? — 6^ is the length of nuclear matter that the dipole traverses 
at impact parameter h for a uniform spherical nucleus of radius i?, and A is 
the mean free path for inelastic dipole-nucleon interactions. Using 



13 



FIG.7 



A = [pa] 



1 



(32) 



with p the nuclear density one has 



(33) 



if we suppose S is purely real. Detailed calculation gives |T0[] 



(34) 



Problem 6(H*): Check (34) for scattering of a dipole on a bare quark. 

Thus we know how to calculate the last three terms in [ ] in (29). What 
about the S'^S term? Graphically this term is illustrated in Fig. 8 for some 
typical elastic and inelastic interactions of the dipoles with the nucleons in 
the nucleus. Let us focus on the final interaction, the one closest to the cut 
(the vertical line) in Fig. 8. One can check that interactions with the x = 
line cancel between interactions in the amplitude and the complex conjugate 
amplitude. 



Problem 7(H): Verify, using S'^S = 1, that the three interactions shown 
in Fig. 9 cancel. You may assume that T is purely imaginary {S = 1 — iT) 
although this is not necessary for the result. 
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FIG. 10 
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FIG. lib. 

Problem 8(H): Again, using S'^S = 1, show that the two interactions shown 
in Fig. 10 cancel. Assume T is purely imaginary. 



The results of problems 7 and 8 establish that one need not consider inter- 
actions with the x = line in evaluating the S'^{x2)S{xi) term in (29). Thus, 
we are left only with interactions on the Xi— line in the amplitude andwith 
the X2-line in the complex conjugate amplitude. The result of problem 8 
allows one to transfer the interactions with the X2-line in the complex conju- 
gate amplitude to interactions with a line placed at X2 in the amplitude 



although this does require that 5* be real. This is illustrated in the equality 
between the terms in Fig. 11a and Fig. lib and reads 
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S''ix2)S{xi) = S{xi — X2). (35) 

Eq.(35) is a beautiful result, but unfortunately it likely has a limited validity. 
It appears to require a maximum of two gluon lines exchanged between any 
given nucleon and the dipole lines as well as the requirement of a purely real 
S-matrix for elastic scattering. Thus using (29), (33) and (35) we arrive atp 

[1 + (.-{2ii-2i2fQl/^ _ e-;^?0|/4 _ g-a;iQ|/4j ^3g^ 



with 



and 



Ql = ^Ql (37) 



Ql = v/^^ P^G, (38) 

where Q5 is the quark saturation momentum and is the gluon saturation 
momentum. In (38) xG is the gluon distribution in a nucleon. Eq.(36) is a 
complete solution to the sea quark distribution of a nucleus in the one-quark- 
loop approximation. 

Problem 9(M-H): Use (30) and (36) to show thatp 

dxjqf + qf)A ^ 

when f /Ql « 1. 

Eq.(39) gives meaning to the idea of saturation as a maximum occupation 
number for, in this case, quarks. As the nucleus gets larger and larger, that 
is as R grows, grows and so there are more and more sea quarks in 
the nuclear wavefunction, nevertheless, the 2-dimensional occupation number 
hits a constant upper bound for momenta below the saturation momentum. 

Problem lO(M-H): Show that 1/2 of (39) comes from the first term in 
[ ] in (36) and that 1/2 comes from the second term with the third and 
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fourth term being small. The second term can be viewed as due to inelastic 
reactions while the first term is the elastic shadow of these inelastic reactions 
of a dipole passing over the nucleus. 

9 Gluon Saturation in a Large Nucleus; the 
McLerran Venugopalan Model^^ 

In order to directly probe gluon densities it is useful to introduce the "cur- 
rent" 

J = (40) 

We shall then calculate the process j + A —>■ gluon {£)+ anything. I shall 
interpret the process in a slightly modified Bjorken frame, and in light-cone 
gauge, while we shall do the calculation in a covariant gauge and in the rest 
system of the nucleus 0. 

For the interpretation we take the momenta of a nucleon in the nucleus 
and the current to be 

P = (P+ t;— >0,0,p) 
2p 

q = (0,0,0, -22;p). 

For small x p is large since P = ^ and this frame is much like an infinite 
momentum frame. If we choose an appropriate light-cone gauge, one that 
eliminates final state interactions, then the transverse momentum and the x 
distribution of gluons in the nuclear wavefunction is the same as the distri- 
bution of produced gluon jets labeled by £ in Fig. 12. 



In order to do the calculation of the spectrum of produced gluon jets we 
carry out a multiple scattering calculation, a term of which is illustrated in 
Fig. 13 H. The calculat ion is simplest to do in covariant gauge. Then 
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FIG. 13 



where xG is the gluon distribution for a single nucleon, and we make use of 
the coordinate space interpretation where it corresponds to the gluon in the 
complex conjugate amplitude being separated from that in the amplitude 
by Ax_L ~ 1/Q- We note that the momentum space distribution of gluon 
produced off a single nucleon, the unintegrated gluon distribution, is given 
by 



dxG 



4^ 



xG{x, 1/x^ 



-il-x 



(42) 



as can be easily verified by integrating the lefthand side of (42) over d?lQ{Q'^ — 
^). Carrying out the 2;— integration in (4) one gets|T3| 
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c?^x e -- (43) 



where Q\ is given in (38). We also note that the e~i^^'sl'^ factor in (41) is 
just 5''''(6)5'(6 + x) = S{x) for a gluon dipole to pass over a length z of nuclear 
material. Thus the derivation of (41) closely resembles that leading to (35) 
but now for gluons rather than quarks. Finally, when |^/<5| << 1 we find 
from (43) 

while for f /g| » 1 Eq.(43) gives the nuclear gluon distribution as simply a 
factor A times the nucleon gluon distribution. Eq.(44) shows that saturation 
is somewhat more complicated for gluons. The factor of A'^ — 1 counting 
the number of species of gluons is expected as is the factor of aN^ in the 
denominator. What is a little surprising is the log factor for which there is 
not yet a good intuitive understanding. Whether this log is a general factor 
or a peculiarity of the present model is not known for sureR T\\. 



10 The Golec-Biernat WiisthofF ModeP^ 

We turn for a while to some phenomenology to see whether there is evidence 
for saturation of parton densities in deep inelastic lepton-proton scattering. 
So far the best way that has been found to approach this problem is through 
a simple model of deep inelastic and diffractive scattering inspired by the 
idea of saturation. We can motivate this discussion by going back to (36) 
and, supposing that such a picture might apply to a proton as well as a large 
nucleus, summing over / and integrating over I and h obtain 



i^2(x,Q^) = / Sx ['dzJ2\i,Ux,z,Q)\' [ d%[l~e-^'^i/% (45) 

Now our Q| naturally depends on the impact parameter b, as indicated in 
(37) and (38) for a nuclear target, however, as an approximation we suppose 
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where Rq will now be taken to depend only on x. Thus 



(46) 



F2ix, Q') = / d\ f rf^E IV^^Afe Q)?{1 - e-^'/<)ao (47) 



which is the formula used by Golec-Biernat and Wusthoff|[T5[| . In addition it 



is then natural to take the diffractive cross section to be given by the shadow 
of the inelastic collisions in which case one replaces cro(l — e- /^^o) = crg(l — S") 
by |(To(l - Sf giving 



F-(x,Q2) = / d-'x f dz z,Q)\'laoil- e-^'/'^of. 

(48) 

Eq.(47) represents a total cross section, (48) represents an "elastic" cross sec- 
tion while the inelastic contribution would have a factor of |o"o(l ~ e^- /^^o) 
replacing the last factors in (47) and (48). 

Golec-Biernat and Wiisthoff include a quark-antiquark-gluon scattering 
term in addition to the quark-antiquark dipole term given by (48) so that 
larger mass diffractive states can also be described. For our purposes this 
is a detail which in any case introduces no new parameters. The model has 
three parameters 

ao = 23mb, Rq^ =< Q| >= {^)^GeV^ (49) 

where A = 0.3 and Xq = 3xl0~^. With these three parameters a good fit to 
low and moderate and low x F2 and F2 data is obtained. In fact the 
fit is surprisingly good over a range of which is remarkably large given 
that there is no QCD evolution present in the model. We shall have to wait 
for further tests and refinements to be sure that the fits are meaningful, but 
we may have the first bit of evidence for saturation effects. The fact that 
< Ql > is in the IGeV^ region is reasonable. 
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11 Measuring Dipole Cross Sections 



(50) 



Refer back to (47). It is easy to check that 

1/Q2 

when Q'^/Q'^g » 1. Thus, although we may view F2 as being given by 
a dipole cross section the size of the dipole is not well determined by 
but, rather, varies between 1/Q and 1/Qs- Thus currently deep inelastic 
structure functions are not well suited for determining dipole cross sections. 
The situation should improve considerably when the longitudinal structure 
function is measured, for in this case the dipole size will be fixed to be of size 
1/Q and will give a direct measure of the dipole cross section. 

Problem 11(E): Use 30 and (45) to derive (50). What is c? 

At present the best place to measure dipole cross sections, and thus to 
see how close or how far one is from finding unitarity limits, appears to be 
in the production of longitudinally polarized vector mesons [p!^ . 



The cross section for 72+ proton — > pl+ proton can be written as 
da^L^P 1 , /• ,0 .0, 



dt 



— I Jd^xd^bJ^ dz^;{x,z){l- S{x,b))tP^4x,z,Qy^-^f (51) 

where t = — is the momentum transfer. The process is pictured in Fig. 14 
and can be viewed in three steps, (i) The virtual photon breaks up into a 
quark- ant iquark dipole of size x with z being the longitudinal momentum 
fraction of the 7* carried by either the quark or antiquark.(ii) The dipole 
scatters elastically on the proton with scattering amplitude 1 — S{x, b) where 
X is the dipole size and b the impact parameter of the scattering, (iii) The 
quark- ant iquark pair then become a p long after they have passed the proton. 
This is the sequence of steps given in (51) where one also integrates over all 
possible dipole sizes and where the integration over impact parameters gives 
a definite transverse momentum. 



Now suppose S is purely real, and define N{Q) = (ipp, ip^*). Then one can 
take the square root of both ides of (51), and after taking the inverse Fourier 
transform one finds 
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FIG. 14 




< S{x,ro,b) >= 1 - J £/^e~^^--^\l^. (52) 

In (52) < S >= {ipp, Sip^) while x denotes the Bjorken-x value of the scat- 
tering and ro(Q) denotes the typical value of the dipole size contributing to 
(51). For (52) to have meaning it is important that the range of dipole sizes 
contributing to (51) not be too large, and this appears to be the case for 
longitudinal p production. This analysis is modeled on the classic analysis of 
Amaldi and Shubert||17| for proton-proton elastic scattering. Thus one can 



estimate the S"— matrix for a dipole of size tq, determined by Q, scattering 
on a proton at a given impact parameter if the data are good enough to 
carry out the integral in (52). The analysis is not model independent as 



one needs to take a wavefunction for the |T^, ^in order to evaluate N. 
However, does not appear to be very sensitive to the choice of wavefunc- 
tion. Also the data are not good enough at larger values of t to accurately 
determine 5* below b ^ 0.3. fm. One finds, for example, at = 3.5GeV 
where tq ^ 1/5/m and for x ^ 10~^ that: (i) S{b ^ 0) ^ 0.5 — 0.7; (ii) 
the probability of an inelastic collision = 1 — S^{b) is considerable at small 
values of b indicating a reasonable amount of blackness at central impact 
parameters; (iii) o'1^i~^^"^°^ ~ 10m6; (iv) is consistent with that found 
in the Golec-Biernat Wiisthoff model. This adds, perhaps, another piece of 
evidence that saturation is approached in the HERA regime for moderate 
values of Q^. 
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12 A Simple Equation for Unitarity; 
the Kovchegov equation^^ 



Consider a (not too high momentum) dipole scattering on a high-energy 
hadron. We suppose the quark- antiquark dipole is left moving while the 
hadron is right moving. Further we suppose that the rapidity, y, of the dipole 
is such that ay << 1 so that one need not consider radiative corrections to 
the wavefunction of the dipole to evaluate the scattering amplitude. We wish 
to study the dependence of the elastic scattering amphtude as one changes 
the relative rapidity of the dipole and the hadron by an amount dY when the 
relative rapidity is Y. Clearly one can view the change dY either as increasing 
the momentum of the hadron and thus allowing its wavefunction to evolve 
further or as increasing the momentum of the dipole. The latter is easier 
to deal with since the dipole is a simple object. When the rapidity of the 
dipole is increased there is a small probability, proportional to dY, that the 
dipole emits a gluon before it scatters off the hadron. We now calculate the 
probability for producing this quark-antiquark-gluon state. 

Since a gluon is emitted either off the quark or off the antiquark we have 
already done the basic emission amplitude, and it is given in (8). It will 
be convenient to work in a basis where transverse coordinate are used rather 
than transverse momenta so one must take the Fourier transform of (8). Thus 
the amplitude for a quark having transverse coordinate Xi to emit a gluon 
having transverse coordinate z, longitudinal momentum k- (Recall that in 
Sec. 2 we were dealing with right movers while here we are concerned with 
left movers.) and polarization A is 



with ipxi^) given by (8) with the replacement kj^ ^ k-. We are now using 
— light-cone gauge to evaluate the left moving quark-antiquark-gluon 
state. Using 




(53) 




(54) 



one gets 
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FIG. 15 



To calculate the probability of a quark-antiquark-gluon state one adds the 
graphs in Fig. 15 to get[E3] 



, 2 2'^ (27r)32A;_^ (z-x^yiz-x^y {z-x^^ (z-x^y' 

(56) 

or 

dn- 2n^'^''^^{x,-znx,-zy ^^^^ 

where we have set dk^/k- = dY and Ci;' = ^ in the large Nc limit where 
the Kovchegov equation will be valid. 

Then the S"— matrix for the quark-antiquark-gluon state to elastically 
scatter on the hadron multiplied by dPr in (57) gives the change in the 
S*— matrix, dS, for dipole-hadron scattering. The result is the Kovchegov 
eg nation ||2l[] 



dS{xi - X2,Y) aNc f (^1-^2)^ 



dY 27r2 7 {xi-zy{x2-zf 

[S{x^ - z, Y)S{z - X2, Y) - S{x^ - X2, Y)l (58) 
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FIG. 16 



and it is illustrated in Fig. 16. We have assumed that the scattering of the 
two dipoles, the quark- (antiquark part of the gluon) and the (quark part 
of the gluon)-antiquark dipoles, factorize when scattering off the hadron. 
This was clear in the model Kovchegov considered where the hadron was 
a large nucleus. This factorization is less obvious in the general case and 
the Kovchegov equation may be a sort of mean field approximation to a 
more complete equation. Also the final term on the right-hand side of (58), 
corresponding to the last two graphs in Fig. 16, give the virtual contributions 
necessary to normalize the wavef unction [^] . The necessity of this last term 



can be seen by considering the weak interaction limit where S ^ 1. Then 
the final term on the righthand side of (58) is necessary to get ^ = when 
^= 1. 



There are two interesting limits to (58). First suppose that S is near 1 
and write S = 1 — iT. One easily finds, keeping only linear terms in T, 

dT{x^-X2,Y) aNc f ^ (^1-^2)^ rr^^ v\ ^t-^ v\] 
ay 11'^ J — zY[x2 — zY 2 

(59) 

which is the dipole form of the BFKL equation. In this case the factorized 
form of the scattering is justified by the large Nc limit and the weak coupling 
approximation. 

The other interesting limit is where 5* is small in which case one need 
only keep the second term on the righthand side of (58) giving 

dS{xi — 2L2^'^) C(^c f (Pz{xi — X2^'^ 



L % ,, g(xi-x,y). (60) 
J uci — zYixo — zY 



dY 2772 J {x^- zY{x2- zf 
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Of course (60) as written cannot be valid. The assumption that S be small 
can be true only when the dipole size is large compared to 1/Qs- Thus we 
should restrict the integration in (60) to the region (x^ — ^.2)^ >> l/Ql as 
well as to the region [xi — z_y, — z)"^ » VQl ^^at the nonlinear term 
in (58) not cancel the linear term. In the logarithmic regions of integration 
one can rewrite the integral (60) as 



dS{x^ - X2, Y) aNc /-fe-a)' 7rd(xi - z) 



dY 27r2 Ji/Q 



giving 



/ , z^S{xi- X2,Y} 61 



dS{x,Y) ^^C(,_if^2 2 



dY TT 

whose solution is 



tn{Q%x')S{x,Y) (62) 



S{x,Y) = e-^^-o'^ '""^^^^)^'^5(x,yo). (63) 
If Ql is exponentially behaved 

QUy) = e^'^^y-'-^^QliY,) 

then 

5(x, Y) = e-f (^)'(^-^«)'5(a;, Fo) (64) 
where Yq should be chosen to satisfy 

QUyo)x^ = 1. (65) 

Eq.(64) is the result found in Ref.23 and argued heuristically already some 
time ago ||24|| , although without an evaluation of the coefficient of the Y'^ term 
in the exponent. 

The Kovchegov equation is an interesting equation for studying scattering 
when one is near unitarity limits. In the next section we shall use a procedure 
very similar to what we have done here to derive an equation whose content is 
presumably equivalent to the Balitsky equation^], a somewhat more general 
form than the Kovchegov equation. The advantage of (58) is its simplicity 
and its likely qualitative correctness in QCD. 
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13 A Simple Derivation of the 
JIMWLK equationt^^ 25-27] 

Over the past seven years or so there has been an ambitious program ded- 
icated to finding appropriate equations for deahng with high density wave- 
functions in QCD. This program has been quite successful and a renormal- 
ization group equation in the form of a functional Fokker-Planck equation 
for the wavefunction of a high-energy hadron has been given by the authors 
of Refs.l4 and 25-27. (JIMWLK). The most complete derivation is given in 
Ref.27 where the equation is written in terms of a covariant gauge potential, 
a, coming from light-cone gauge quanta in a high-energy hadron. Here we 
give an alternative simple derivation ||28||. 

We can imitate this mixture of gauges used in Ref.27 by taking Coulomb 



gauge p9| which has a gluon propagator 



Dapik) = -^[9af^ ^^-^ ] (66) 

where N ■ v = Vq for any vector v. Suppose the propagator has k'^ » 
k^, fc^ and connects two highly right moving lines. Then the Coulomb gauge 
propagator is equivalent to the = light-cone gauge propagator 

Dap{k) = -^[gap ] (67) 

Problem 12(E): Show that D , and D^i as given in (66) and (67) agree 

when k^ >> k^,k'^. For a right moving system these are the important 
components of the gluon propagator. 

Similarly for a left moving system Coulomb gauge is equivalent to A_ = 
gauge while for gluon lines which connect left moving systems to right moving 
systems the dominant component in (66) is = l/(fc)^ which looks like 
covariant gauge when fc^, /c^ « k^. 

We are going to consider the scattering of a set of left moving quanta, 
quarks and gluons, on some high-energy, right moving hadron. These quarks 
and gluons may be parts of a hadronic wavefunction which are frozen in the 
passage over the right moving hadron or they may come from a current as 
in our discussion of deep inelastic scattering given in Sec. 8. For simplicity 
we shall limit our discussion to left moving quark and antiquark lines, but 
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this is simply to avoid too cumbersome notation. Then a left moving quark 
interacting with the right moving hadron can be represented by 

(x) — P exp{ig f dx-A+{x,X-) (68) 

J—oo 

where we have taken a;+ = and fixed the left moving quark to have trans- 
verse coordinate x. Except for the change of right moving quark to left moving 
quark (68) is the same as the operator in the matrix element in (26) where 
we showed how quarks could be identified with Wilson lines in the fundamen- 
tal representation. By taking gauge invariant combinations of V^'s and l^^'s 
we can form observables which depend on and which correspond to the 
scattering of quite general left moving systems on the right moving hadron. 
We denote a general such observable by 0(A+). 

Although we have put the integration in (68) exactly on the light-cone 
we in fact are going to assume that the left moving observable has rapidity y 
obeying ay < < 1 so that transverse gluons are unlikely to be emitted by the 
left movers allowing us to identify the left moving system at x_ = ±oo. If 
the relative rapidity of the scattering is Y then we imagine that a(Y — y) pa 
aV >> 1 so that the right moving hadron has, in general, a wavefunction 
including many gluons. If the right moving hadron has momentum p then 
the scattering amplitude is 

< O >y=< p\0\p >= j D[a{x,x^)]0{a)WY[a\ (69) 
where the weight function Wy is given by 

Wvia] = J D[A^]S{A+ - a)S{F{A))AF[A]e'^^^^ (70) 

where D indicates a functional integral, F is a gauge fixing, and A^j^ is the 
corresponding Fadeev-Popov determinant times an operator which projects 
out the state \p > initially and finally. We suppose Wy is normalized to 

J D[a]Wy[a] = 1. (71) 

Now consider the F— dependence of < O >y . From (69) one can clearly 
write 

^<0>y^J D[a]0{a)^Wy[a]. (72) 
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(a) (b) (c) 

FIG. 17 

However, one can equally well imagine calculating <0 >y by evaluating 
the change of the left moving system, due to an additional gluon emission, 
as one increases the rapidity of the left moving system by an amount dY. 
The change in the left moving state is given by the gluon emissions and 
absorptions shown in Figs. 17 and 18 where spectator quark and antiquark 
lines are not shown. In Fig. 17 we have assumed that a gluon connects a 
quark and an antiquark line. This is for definiteness. We equally well could 
have assumed a connection to two quark lines or to two antiquark lines. 
The vertical line in the figures represents the "time," x_ = 0, at which 
the left moving system passes the right moving hadron. This view of the 
F— dependence of < O >y in terms of a change of rapidity of the (rather 
simple) left movers is in the spirit of work previously done by Balitsky[5], 



KovchegovpHI, and Weigertp6|. 



We shall examine in some detail the graphs of Fig. 17 before stating the 
complete result including the graphs of Fig. 18. Begin with the graph shown 
in Fig. 17c. We do the calculation in A- = light-cone gauge, which for left 
movers is equivalent to our Coulomb gauge choice. This graph is exactly the 
same as has been calculated in Sec. 11 except for the fact that the quark and 
antiquark lines are not necessarily in a color singlet. The result is the same 
as I the first term on the righthand side of (56) except for the color factors 
which we put in separately. The result is 
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(a) 



(b) 



(c) 



FIG. 18 



,2^ (^1 - Z) ■ {X2 - Z) 



(73) 



in going from O to The additional factors, as compared to the corre- 
sponding term in (56), are T*^ which comes to the right of V^{xi) because 
the emission off the Xi— line is at early values of x_, the T'' which comes to 
the right of V{x2) because the absorption on the a;2— line is at late values of 
and the factor Vcd{z.) giving the interaction of the gluon with the hadron 
as a Wilson line in the adjoint representation. 
Now write 



late values of x_. Indeed, we may view the graph in Fig. 17c as being given by 
the "mnemonic" graph shown in Fig. 19 where the adjoint line integral starts 
at large positive values of a;_ and proceeds to large negative values of a;_ at 
a transverse coordinate Xi then back again to large positive values of x_ at 
a transverse coordinate z. The result of the graph of Fig. 17c then is 



Now it is straightforward to add in the other terms of Fig. 17 and those 
of Fig. 18 to get 





■^{-V\x,)V{z)U 



(75) 
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X _ 



b 



FIG. 19 



V\x,)(^V{x,) ^ asil I d'xd'yr^^l—- -+ / d'xu^ ^ 



2 J -^6a''{x,x^)6a\y_,y_) J - (5a" (x , x„ ) 

■V\x{)® V{x^) (76) 

where, if one takes large positive values of x_ and the functional deriva- 
tives in (76) simply insert color matrices in the appropriate places, as, for 
example, in (75). rj and v are given by 



9'€y = 4 / 



d?z {x — z)-{y — z) 



An'^ {x — zy {y — z)'^ 
■ {1 + V^{x)Viy) - V\x)Viz) - V\z)Viy)}ab (77) 



and 

£z 



9< = 2^J j^^^tr[T^VHx,Mz)]. (78) 

For a general scattering one simply replaces ^^(x^) ® V{x2) in (76) by 
0{a). After multiplying by PVy[Q;] and integrating over D[a] one gets|| 
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where an integration by parts in a has been done on the righthand side of 
(79). To the extent that the 0{a) form a complete set of observables, and it 
is not clear how close this is to being true, one can equate the integrands of 
(79) and arrive at the JIMWLK equation 



The exact values of x- and y- appear to have some arbitrariness as dis- 
cussed in some detail in Ref.28. Eq.(80) is an elegant equation of a functional 
Fokker- Planck type the nature of whose solutions is now under investigation. 
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